Role of electrostatic interactions in the assembly of empty spherical viral capsids 
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We examine the role of electrostatic interactions in the assembly of empty spherical viral cap- 
sids. The charges on the protein subunits that make the viral capsid mutually interact and are 
expected to yield electrostatic repulsion acting against the assembly of capsids. Thus, attractive 
protein-protein interactions of non-electrostatic origin must act to enable the capsid formation. We 
investigate whether the interplay of repulsive electrostatic and attractive interactions between the 
protein subunits can result in the formation of spherical viral capsids of a preferred radius. For this 
to be the case, we find that the attractive interactions must depend on the angle between the neigh- 
boring protein subunits (i.e. on the mean curvature of the viral capsid) so that a particular angle(s) 
is (are) preferred energywise. Our results for the electrostatic contributions to energetics of viral 
capsids nicely correlate with recent experimental determinations of the energetics of protein-protein 
contacts in Hepatitis B virus [P. Ceres and A. Zlotnick, Biochemistry 41, 11525 (2002).]. 
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I. INTRODUCTION 

Many types of viruses can spontaneously assemble 
from the proteins that make the viral coating (capsid) 
and the viral genetic material. This was first demon- 
strated in the work of Fraenkel-Conrat and Williams 1 
where they showed that fully infectious tobacco mosaic 
viruses can be spontaneously reconstituted from the indi- 
vidual proteins that make the viral coating and the viral 
RNA molecules. 

It may be argued that the precision of the viral as- 
sembly is guided by an interplay between the proper- 
ties of the RNA/DNA (bio-polyelectrolyte) and the cap- 
sid proteins, such as the total bio-polyelectrolyte length, 
the effective charge on the proteins and on the bio- 
polyelectrolyte, and possibly some specific interaction 
acting between the RNA/DNA and the proteins. How- 
ever, in some types of viruses, the genetic material is 
not necessary for the capsid assembly and empty viral 
capsids can be assembled in the absence of viral bio- 
polyelectrolyte^i 3 .'^^, at least when the amount of salt 
in the solution is large enough 6 . It is rather intrigu- 
ing that the thus assembled capsids are often highly 
monodisperse^i 3 . which suggests that there is some reg- 
ulating factor, independent of the viral RNA/DNA, that 
favors capsids of particular size. This is the effect that 
we investigate in this article and the motivation for our 
study. 

In principle, viral proteins can assemble in a variety of 
capsids with different shape and size 5 -^. We shall con- 
centrate on the nearly spherical viral capsids whose struc- 
ture can be described within the so-called Caspar-Klug 
quasiequivalent construction-. In more mathematical 
terms, the viral capsids that we consider are icosadeltahe- 
dral, i.e. they can be mapped onto triangulations of the 
sphe re with the icosahedral "backbone" (see e.g. Refs. 
l8l ll0lfTll for more details). The different triangulations 



can be described using the notion of the so-called T- 
number. The number of protein subunits (N, which may 
also be protein dimers, such is e.g. the case in Hepatitis 
B virus) in a quasiequivalent viral capsid is 

N = 60T, (1) 

where T = 1,3,4, 7,9, i.e. T = h 2 + hk + k 2 and h 
and k are non-negative integers^. 

In view of the potential polymorphism of the viral pro- 
tein assemblies which is also observed experimentally, es- 
pecially in some viruse o 6 ' 12 : 13 , it is rather surprising that 
empty capsids of most simple viruses precisely assem- 
ble in capsids whose T-number is the same as in fully 
functional capsids that contain the viral RNA/DNA 2 ^ 3 -. 
Even when several differently sized empty capsids do 
form^ii 3 -, these represent only a tiny subset of all of 
the imaginable capsid formations. Since the viral pro- 
teins have certain charge in the solutio n 15 : 16 , it is tempt- 
ing to assume that the repulsive electrostatic interac- 
tions compete with some other attractive interactions 
between the viral proteins in such a way that the to- 
tal free energy of the capsid is minimal exactly for the 
capsid of the observed radius. The attractive interac- 
tions between the viral proteins could be of different ori- 
gins (e.g. van der Waals interactions, hydrophobic in- 
teractions or chemical bonding). However, it has been 
experimentally demonstrated that the binding energy 
of the two viral proteins increases with temperature 3 , 
which strongly suggests that the attractive protein inter- 
actions are dominated by either the zero-frequency term 
of van der Waals interactions 1 ^ and/or by hydrophobic 
interactions^. They should thus be rather local and pro- 
portional to the area "buried" in the protein-protein con- 
tacts. 

The competition between hydrophobic and electro- 
static interactions in viral capsids has been theoretically 
investigated befor o 15 : 18 ' 20 . The emphasis there, however, 
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was mostly on kinetics of viral assembly, i.e. on the clas- 
sical nucleation theory and the mass action law. Further- 
more, the electrostatic interactions were modeled as the 
asymptotic form of the Debye-Hiickel (DH) approxima- 
tion, in the regime where R/Xdh 3> 1, where R is the 
capsid radius and Xdh the DH screening length that is 
proportional to inverse square root of the salt concentra- 
tion. It is the aim of this article to reexamine in more 
detail the role of electrostatic interactions in the capsid 
assembly, to investigate whether they can act to regulate 
the capsid size (radius) and whether they can prevent the 
capsid formation. To this end, we shall introduce a model 
of a viral capsid that is more realistic and more complex 
than the ones used in previous studies. We shall go be- 
yond the DH approximation and consider all regimes in 
R/Xdh- We shall also estimate the contribution of the 
electrostatic interactions to the total binding energy of 
the capsid. 

There are two different models of capsid electrostat- 
ics that we consider. Section [H] describes the simpler of 
these models that treats the capsid as uniformly charged, 
infinitely thin sphere whose charge density is fixed by the 
total charge on the capsid protein in the solution. This, 
rather simplified, model of capsid allows us to relate our 
results with those previously published 15 . The aim of 
Section [II] is twofold. First, we shall clearly demonstrate 
that the previously derived expressions for the capsid 
energetics^ have limited validity, even in physiological 
conditions. The theory and results presented in Section 
Hll also constitute a good prelude for a more elaborate 
electrostatic model of a viral capsid presented in Section 
IIIII There, the capsid is modeled as a dielectric medium 
contained in between the two infinitely thin spheres, each 
of which has a certain prescribed charge density. In Sec- 
tion IIVI we attempt to relate our theoretical predictions 
to experimental results on the viral energetics^. Section 
[V] discusses limitations of our models of viral capsids. 



II. VIRAL CAPSID AS AN INFINITELY THIN, 
UNIFORMLY CHARGED SPHERICAL SHELL: 
MODEL I 

A. Free energy of empty capsids in salty solution 

We first approximate the icosadeltahedral capsid as the 
perfect sphere of radius R whose charge is uniformly dis- 
tributed on the surface, so that the surface charge density 
is a. The theory of elasticity applied to icosadeltahe- 
dral shells predicts that the viral capsids are aspherical, 
the more so the larger their mean radius (this is also in 
agreement with experimental data) 11 ' 21 . The asphericity 
arises from the so-called "buckling" of the capsid around 
the pentameric protein aggregate s 11 ' 21 . The experimen- 
tal data on virus shapes 2 ^ suggests that the virus surface 
may in fact be very corrugated, but it is difficult to as- 
sess the corrugation and asphericity of the corresponding 
protein charge distribution. For our purposes we neglect 



the capsid asphericity which is, on the basis of continuum 
theory of Lidmar et ai^I, expected to be small especially 
in viruses of small radii (the validity of the theory may, 
however, be questionable for small viruses; see Sec. IVf . 
Under this approximation, the charges on the capsid, to- 
gether with the (monovalent) salt ions in the solution 
whose bulk concentration is Co, give rise to spherically 
symmetric electrostatic potential &(r), so that the prob- 
lem is effectively one-dimensional. 

The free energy of the system (proteins and salt ions 
in the solution) in the mean-field (Poisson-Boltzmann) 
approximation can be expressed as a functional of the 
electrostatic potential as 2 ^ 

F = J d 3 r [f e i(r) + f ton s{r)} + F boundary , (2) 

where 

fa(r) = ec +$- ec -$-^|V$| 2 , (3) 

and 

fions{r) = £{±[c i (r)]n C i (r)-c i (r) 

(4ln4-4)] _ M *[ c *(r)-4]}, (4) 

where Fb oun dary is the boundary contributions arising 
from the discontinuity of the potential at the capsid 
[see Eqs. ([7]) and © below]. The free energy is the 
sum of the electrostatic energy of charge in the poten- 
tial (fei), and the salt ions configurational/cntropy con- 
tribution (/ions). This ansatz is correct as long as the 
counterions are of low valency and/or the charge den- 
sity on the capsid is not too large, leading to the s.c. 
weak coupling-regime which is properly captured by the 
Poisson-Boltzmann theory^. 

In the equations above, e is the electron charge, e r is 
the relative permittivity of the solvent (water, e r = 80), 
eo is the vacuum permittivity, f3 — {ksT)~ l , T is the tem- 
perature (T=300 K), and ks the Boltzmann constant, 
c + (r) and c~ (r) are the concentrations of + and — ions 
in the solution and /i + and fi~ are their respective chem- 
ical potentials. We shall consider dilute capsid protein 
solutions and shall thus neglect the influence of counte- 
rions on the potential. This also means that the bulk 
concentrations of positive and negative ions are equal, 

4 = c o = c o- 

The variation of the functional Eq. ([2]) with respect to 
fields $ and c l yields the following equations: 

c^r) = c exp[T-e/3$(r)], (5) 

and 

e erV 2 $(r) = 2ec sinh[/3e$(r)]. (6) 

Equation (|6|) is the Poisson-Bolzmann equation for the 
potential field. Its mean field character is a consequence 
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of approximating the salt ions as the ideal gas [Eq. Q]. 
One additionally needs to specify the boundary condition 
of the <!> field. For the assembly of virus-like particles, the 
prescribed surface charge density (a) is the most realis- 
tic boundary condition and regardless of the capsid size, 
i.e. the number of protein subunits it contains, the sur- 
face charge density of capsid should be the same. We 
shall also assume that the charge on the particular pro- 
tein does not depend on the salt concentration (but see 
Ref. [H for the dependence of protein charge on the pH 
value of solution) and will thus exclude charge regulation 
boundary condition from our considerations 17 . Requiring 
that the surface charge density of the bounding sphere is 
fixed, yields 



<9$(r) 



dr 



d$(r) 



r=R- 



dr 



= 7T- ^ 

r=R+ e 0tr 



The derivative of the potential (or electric field) displays 
a discontinuity at the bounding sphere (see Fig. [1]). This 
is of course a consequence of the fact that there is a charge 
(Q = 4naR 2 ) on the sphere. The electrostatic energy of 
this charge is equal to Fb oun( iary = Q$(R) and the total 
electrostatic free energy of the system is thus 



F = / [fel(r) + nons(r)]d z r + Q$(R), 



(8) 



Note that in the limit Co — > (the Coulomb limit), the 
functional reduces to 



F 



2 



J \V$\ 2 d 3 r + Q${R) 



(9) 



which is simply the self-energy of the (completely un- 
screened) charge on the sphere. 



B. Numerical solutions of the model and analytic 
approximations in the Debye-Hiickel limit 

Equation ^ is a nonlinear differential equation which 
we solve numerically. We first discretize the radial coor- 
dinate both within the capsid and outside it, so that the 
total number of points is typically around 600 25 . The 
intervals in radial coordinates are not the same in the 
two regions, and have to be chosen so that the potential 
at the point with largest radial coordinate (outside the 
sphere) is very nearly close to zero. At each of the points 
we consider the value of the electrostatic potential. The 
sum of square deviation of the potential from the values 
required by the differential equation (JSJ) at the discrete 
set of radial coordinates is minimized 24 until a desired 
numerical accuracy is achieved. Similar numerical proce- 
dure has also been used in Ref. 0. 

Additional insight can be obtained in the regime when 
e/3|$| <C 1. In that case (the DH approximation) the 
equation can be linearized and analytically solved. The 



solution within the capsid is 

Qsmh(n DH r) 



$(r) = 



4irreoe r KDHR [sinh(KDjji?) + cosh(K£>#i?)] 



and outside the capsid 

Qexp[-K DH (r - R)] 



$(r) 



Airrtoe r KDHR [1 + cot\i{n dh R)] ' 



(10) 



(11) 



where km = 1/Xdh = \J (2e 2 co)/ (eoe r ksT) is the in- 
verse DH screening length. Comparison of the numerical 
and the DH solution from the above two equations is 
shown in Fig. [TJ Note how the DH approximation fails 
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FIG. 1: Numerically exact solution for the potential $(r) 
(symbols) and the DH solution (lines). The (mono-valent) 
salt concentrations are 100 mM and 10 mM and the capsid 
radii are _R=11.07 nm and 7?=7.02 nm in panels (a) and (b), 
respectively. The surface density of capsid charge is a = 0.4 
e / nm 2 . 

when the potential at the capsid acquires high values, 
as is the case in low-salt solutions. The surface charge 
density chosen in the displayed result (a = 0.4 e/nm 2 ) 
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should be realistic for the assembly of virus-like particles 
at neutral pHl£i±&I£. 

The free energy of the assembled capsid in salty solu- 
tion can be calculated from Eq. ([5]) once the electrostatic 
potential has been obtained. It is instructive to exam- 
ine the Coulomb case first. In this situation, the total 
free energy is simply the electrostatic self-energy of the 
capsid, 



F c 



oulomb 



Q$(R) _ 2R 3 wo- 2 



£0 e r 



(12) 



For a given surface charge density, the free energy should 
scale with the third power of the capsid radius. If one 
assumes that the attractive protein-protein interaction 
(whatever its source may be) is local, one can write 



2 ' 



(13) 



where N is given by Eq.|T]), m is the mean number of 
contacts that a particular protein makes with its neigh- 
bors (m ~ 3), and 7 is the negative free energy of the 
protein-protein contact (7 > 0). Thus, the contribution 
of F a dhesion to the total capsid energy should scale as 
— R 2 , since the attractive interaction is simply propor- 
tional to the number of protein subunits. The combi- 
nation of the unscreened electrostatic self-repulsion and 
the attractive adhesion interaction should thus have a 
minimum at some radius. This observation could appar- 
ently explain the origin of the monodispersity of capsid 
radii without the need to introduce e.g. the spontaneous 
curvature contribution to the total energy of the capsid. 
However, as shall be shown below, these considerations 
profoundly change when one considers the role of salt in 
the capsid assembly. 

The free energies of the capsids in the solutions 
with several different mono-valent salt concentrations are 
shown in Fig. O It is obvious that even for quite low salt 
concentrations (cq =1 mM) the Coulomb prediction is 
not valid. Furthermore, even the functional dependence 
of the free energy on the capsid radius is different from 
the prediction of Eq. (TT2"]) . In the DH approximation, 
the free energy of the system can be shown to be given 
by 



F, 



DH 



QHR) 



(14) 



somewhat deceivingly looking like the expression in 
Eq. (|12p since the potential at the capsid, $>(R) is to 
be obtained as the solution of the linearized Poisson- 
Boltzmann equation, i.e. it is profoundly influenced by 
salt ions. Using Eq. (|10p. one obtains that 



F. 



2na 2 R 2 



DH 



KDH£0tr[l + COth(K dhR)] 

When kjjhR 3> 1, the equation reduces to 
hm F DH = , 



(15) 
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FIG. 2: Free energies of the system (capsid + salt) as a func- 
tion of the capsid radius for several different salt concentra- 
tions as indicated in the body of the figure. The capsid sur- 
face charge density is a = 0.4 e/nm 2 . The symbols denote the 
numerically obtained (exact) results, while the dashed lines 
are the DH approximations to the free energy (the length of 
dashes increases with the salt concentration). The full line de- 
notes the free energy in the Coulomb regime, Eq. ((12}. The 
points denoted by A and B correspond to potential profiles 
that are shown in panels (a) and (b) of Fig. [T] respectively. 



which shows that for kdhR 3> 1 the free energy scales 
with the second power of the capsid radius in clear 
contrast with the prediction obtained in the Coulomb 
regime, Eq. (|12l) . Exactly the same relation was ob- 
tained in Ref. [l5| [see their Eq. (4)] by using a differ- 
ent approach, but note that its regime of validity is lim- 
ited by the fact that (i) it was derived by linearizing the 
Poisson-Boltzmann equation, and that (ii) it holds only 
when kdhR ^ 1- Thus, for very low salt concentrations, 
the validity of Eq. (fTS]) is severely limited, but even for 
moderate salt concentration of Co = 10 mM (with a cho- 
sen as before), the DH approximation overestimates the 
free energies at R ~ 10 nm by about 50 %. For larger 
surface densities, the DH approximation is worse, and 
for sufficiently large capsid charge density it becomes er- 
roneous even in the physiological salt regime (cq ~ 100 
mM). This is illustrated in Fig. [3] which displays the free 
energies of the assembled capsids at Co = 100 mM for sev- 
eral different values of the effective capsid charge density 
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(cr = 0.4,0.8,l-2,l-6 e/nm 2 ). 
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FIG. 3: Panel (a): Free energies of the system (capsid + salt) 
for Co =100 mM as a function of the capsid radius for sev- 
eral different capsid charge densities as indicated in the body 
of the figure. The symbols denote the numerically obtained 
(exact) results, while the dashed lines are the DH approxima- 
tions to the free energy (the length of dashes increases with 
the capsid charge density). Panel (b): Comparison of nu- 
merically exact results (symbols) with the DH approximation 
(dashed line) for capsid radius R =20 nm. 

Interestingly enough, in the opposite regime of capsid 
sizes, i.e. when kdhR *C 1, the DH expression for the 
system free energy reduces to 

lim F DH = tt—— — 

k dh R<s.i £o£r(l + kdhK) 



(17) 



which clearly tends toward the (unscreened) Coulomb 
expression for the free energy when kdhR — > 0. It should 
be noted here that there is in fact an upper limit on Xdh 
set by the concentration of counterions that the proteins 
release in the solution, so that the limit Xdh — ► oo (or 
k dh — ► 0) should be considered with caution. 

The numerically exact solutions also display the scal- 
ing laws predicted by the DH approximation, i.e. the 
free energy scales as i? 3 for KdhR *C 1 and as R 2 for 



kdhR ^> 1 (see particularly the cq=1 mM case in Fig. 
[5]). One important difference though, is an onset of R 2 
scaling behavior for smaller values of R than predicted 
by the DH calculation, especially when the salt concen- 
tration is low. The same scaling behaviors have been 
obtained in a DH approximation for a more complicated 
model of a spherical viral capsid whose charge distribu- 
tion is nonuniform and carries a signature of its icosahe- 
dral symmetry 26 . 



III. VIRAL CAPSID AS A DIELECTRIC 
MEDIUM CONTAINED IN BETWEEN THE 
TWO UNIFORMLY CHARGED INFINITELY 
THIN SPHERICAL SHELLS: MODEL II 



In the electrostatic model from the previous section we 
have assumed that the capsid is infinitely thin and uni- 
formly charged. The proteins of viral capsids carry a dis- 
tribution of charge, both of positive and negative sign 27 , 
the negative charges being typically localized on the out- 
side and the positive charges on the inside of the cap- 
sid (this trend is significantly more pronounced in DNA 
viruses 2 ^) . The presented model could thus be viewed as 
a smeared lowest order (monopole) in multipole expan- 
sion of the fields generated by the protein charge distri- 
bution. In order to investigate whether the R 2 behavior 
of the electrostatic free energy of the capsid is a conse- 
quence of the simplicity of model I, in this section we 
adopt a more complex model of the capsid, as illustrated 
in Fig. [U The capsid is approximated by the shell of fi- 
nite thickness, S whose dielectric permittivity is e p . The 
shell is permeable to water and salt ions from the outside, 
but no ions are allowed in the capsid material. The inner 
and outer surfaces of the shell carry the surface charge 
densities of o\ and 02, respectively. The thickness of the 
dielectric layer S in our model should not be identified 
with the capsid thickness, although they are clearly re- 
lated. The typical capsid thicknesses are of the order of 
2 nmS 2 -, and this should be the upper bound for the S. 

It is again instructive to obtain some analytical limits 
for the energetics of viral capsids and to this end we solve 
the problem in the Debye-Hiickel approximation. The an- 
alytic results for the electrostatic potential are somewhat 
cumbersome and are thus summarized in the Appendix. 
The capsid free energy is in the DH approximation given 
by 



F. 



DH 



0$ 

2 



27T [crii? 2 $(i?) 



a 2 {R + S) 2 ^(R + S)] . 

(18) 

The number of parameters in this model is large, but 
since we are interested in the assembly of viral capsids, 
we seek for the expression for the free energy in the limit 
when KdhR 3> 1, 6 -C R, and e r > e p (e p is of the order 
of 5 for uncharged portions of the protein assemblies, de- 
pending also on the type of protei n 28 ! 29 ). In this regime, 
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FIG. 4: An illustration of the electrostatic model of the viral 
capsid. The figure represents the cross-section of the assem- 
bled empty viral capsid with parameters of the model de- 
noted. Salt ions are represented by small spheres. The inte- 
rior of the capsid contains water and salt ions, but the salt 
ions are not present in the viral capsid shell. 



the free energy simplifies to 



limFoH = 2ir ep( - ai 



02 J 



e r (a 2 



(19) 

Note that when (5 = the above equation reduces to Eq. 
(JTHJ) with a — <j\ + o-2- At physiological salt concen- 
trations, kuh ~ 1 nm -1 , so that kdhS ~ 1, assuming 
that S is not much smaller than 1 nm. Assuming that 
e r 3> 2e p , a simpler expression for free energy of viral 
capsids is obtained: 



lim F DH = 



(20) 



Note that the above equation reduces to Eq. (|T5|) when 
a 1 = 02 = c/2, but for any other combination of a\ 
and <72 such that a\ + 02 = er, the free energy is larger. 
Interestingly enough, the limiting form of the DH approx- 
imation for the free energy in Eq. ([20]) does not depend 
on e p as long as e p <C e^S. Additionally, in this limit, 
the free energy does not depend on signs of o\ and 02- 

In any case, in the regime typical for viruses, we ob- 
tain the R 2 behavior of the free energy, as in model I. 
To check the validity of the DH approximation for this 
model, and whether the R 2 dependence of the free en- 
ergy is preserved in the exact solution of the problem, we 
have again numerically solved the full nonlinear Poisson- 
Boltzmann equation. The results of these studies are dis- 
played in Fig. [5l which represents the comparison of the 
DH and numerically exact solutions for the electrostatic 
potential, and in Fig. [6] which displays the electrostatic 
energies of viral capsids. 

It is obvious that the DH approximation may be en- 
tirely inapplicable to obtain reliable estimates for viral 
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FIG. 5: Electrostatic potentials obtained numerically (sym- 
bols) and in the Debye-Hiickel approximation (lines). The 
arrows denote positions of the inner and outer side of the cap- 
sid. The parameters of the calculation are e p = 5, S = 1 nm, 
and R — 20.15 nm. Panel (a): ai = 0.5 e/nm 2 , 02 = —0.1 
e/nm 2 , c = 100 mM. Panel (b): ai = 0.6 e/nm 2 , a 2 = -0.2 
e/nm 2 , c = 10 mM 



energetics, especially in the low to moderate salt regime. 
This is most easily seen by comparing the potentials ob- 
tained numerically and in the DH approximation in the 
low salt regime displayed in Fig. [SJd) (note the scale of 
the potential). Note that the R 2 dependence of the free 
energy, predicted by the DH approximation, is confirmed 
by the numerical results in the limit appropriate for viral 
capsids. Intriguingly, inspection of the data obtained by 
numerically solving the problem reveals that the propor- 
tionality of free energy to square of capsid radius holds 
much better than predicted by the DH results, espe- 
cially in the low-salt regime, and even in regime where 
k-dhR ~ 2, at least for the parameters chosen in the 
presented results. This is similar to what was observed 
in model I. Note that the magnitudes of free energy ob- 
tained in model II for a x +a 2 = 0.4 e/nm 2 (2000 to 7000 
ksT at R = 20 nm, cq = 100 mM, depending on o\ and 
02 ) are comparable but always larger than those obtained 
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FIG. 6: Electrostatic free energies of the system (capsid + 
salt) obtained numerically (symbols) and in the Debye-Hiickel 
approximation (dashed lines) for three different combinations 
of the inner and outer surface charge densities, o\ and 02 as 
denoted in the figure. The parameters of the calculation are 
e p = 5 and 5 = 1 nm. Panel (a): Co = 100 mM. Panel (b): 
Co = 10 mM. Full thick lines represent functions that are pro- 
portional to R 2 . The electrostatic potentials corresponding 
to points denoted by A and B are presented in panels (a) and 
(b) of Fig. [S] respectively. 



in model I for a = 0.4 e/nm 2 (« 1800 k B T at R = 20 
nm, Co = 100 mM), in agreement with Eq. (f2"0"|) and the 
discussion following it. 



IV. APPLICATION OF THE RESULTS TO 
VIRUSES 



A. R dependence of free energy and preferred 
mean curvature of the capsid 

In the regime of radii typical for the viral capsids (i? ~ 
10 nm), R 3> \dh limit is clearly appropriate for the 
physiological salt concentrations (cq ~ 100 mM) and thus 
F oc R 2 . The total free energy, i.e. the sum of free en- 



ergies in Eqs. (|SJ and (fl3|) is also proportional to R 2 , 
which means that the free energy per protein is a con- 
stant. This is confirmed by both of the models we con- 
sidered. Since the effect related to entropic differences 
between the (free) proteins in the solution and those 
assembled in capsids are expected to be secondary for 
(large) protein capsids 3 —, the electrostatic interactions in 
combination with the simple expression for the adhesion 
free energy [Eq. (fl~3|) ] cannot explain the occurrence of 
the preferred capsid radii. This is simply because irre- 
spectively of the capsid size, i.e. its T number [see Eq. 
([I])], the free energy per protein remains constant. Thus, 
there is no particular capsid radius that is preferred en- 
ergywise. This is one of the main results of our study. 

A preferred mean curvature of the capsid (Hq) which 
could stem from the curvature of the protein-protein con- 
tacts would enforce a particular capsid radius - this would 
mean that the (non-electrostatic) "adhesion" energy in 
Eq. (fT5)) could be reformulated as 



Nm 



dhesion 



j + uj[H q - — 



(21) 



at least when \Hq — <C 1. Similar views have been 
recently proposed^ 3 -!, but the original idea goes back to 
Caspar and Klug 7 . Our results clearly support the idea 
of the preferred curvature in viral capsids, a subject that 
has recently received much attention^ ' 11 ' 21 ' 31 . If the at- 
tractive interaction has two or more minima in the spon- 
taneous curvature, arising from several stable protein- 
protein conformations, one can also expect formation of 
several differently sized capsids^. It should also be kept in 
mind that allowed radii of the capsid are in fact discrete, 
since the sphere can be triangulated by protein subunits 
only for certain total number of proteins [see Eq. ([Tjl], 
The allowed capsid radii are thus R = ^15TAq/tt. 

Several articles have recently questioned the problem 
of pressure that acts within the self-assembled RNA 
viru a 18 ' 33 . In this respect loosely packed RNA viruses are 
expected to be very different from the double stranded 
DNA viruses whose capsids are known to withstand very 
high internal pressure s 11 ' 3 ^ 35 . It is important to note 
here that the fact that the empty capsids precisely as- 
semble at radii given by R — Hq , which is the same 
as the radius of the capsids filled with the viral genetic 
material (polyelectrolyte) , means that the pressure that 
acts on the capsid arises solely from the polyelectrolyte 
self interaction and its interaction with the capsid^ 3 -. In 
other words, the repulsive electrostatic self-interactions 
in empty capsids can expected to be exactly counteracted 
by the angle-dependent adhesive interactions in Eq. (j2"T]) 
- the first derivative of the sum of these interactions with 
respect to capsid radius (i.e. pressure) is zero. 
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B. Influence of electrostatic interactions the 
strength of protein-protein contacts 



Electrostatics and the inside/outside 
asymmetry of the viral capsids 



The capsid total free energy is negative only for suffi- 
ciently low values of the capsid charge density, i.e. when 



2 2mjAoeoe r 

<7 < 



A 



DH 



(22) 



where Aq is the mean area of the protein subunit in the 
capsid [this estimate is based on the model I; see Eq. 
(|12|)]. This shows that for sufficiently large Xdh the re- 
pulsive electrostatic interaction dominates the energetics 
of the assembly. In solutions of low salt concentrations 
the capsid assembly should be thus inhibited. For suffi- 
ciently small values of adhesion constant 7, our calcula- 
tions predict that there is a critically low salt concentra- 
tion, Cq t at which the assembly does not take place. As- 
suming that the critical salt concentration is high enough 
so that the limit R/Xdh (cg r ) 3> 1 is still satisfied, one 
obtains that cg r = (27r 2 a 4 R 4 k B T)/(e e r N 2 m 2 -f 2 e 2 ). 

Although the presented electrostatic models cannot ex- 
plain the preferred radii of the empty viral capsid, they 
can give some clue with regard to the recently observed 
strengthening of the protein-protein contacts in empty 
Hepatitis B capsids in high salt concentrations^. Namely, 
at high salt, the repulsive protein-protein electrostatic in- 
teractions are weakened as is clearly demonstrated in Fig. 
[2j Ceres and Zlotnick observed that the contact energy 
[7, see Eq. (fl"3]) ] between the capsid proteins of Hepatitis 
B increases from -5.3 k B T at co=150 mM to -6.9 fc^T at 
co=700 mM^. The total increase in energy of the cap- 
sid is thus 240(-6.9 + 5.3) k B T = 384 k B T, since the 
capsid is made of 120 copies of a tetravalent (m = 4) 
protein dimer Cpl49^ [see Eq. (fi"3])]. Our calculations 
based on model I predict that for er=0.4 e/nm 2 , the cap- 
sid free energies at i?=15.6 nm (which is the mean radius 
of Hepatitis B capsid^), calculated from Eq. J2]) are 850 
k B T and 415 k B T for c = 150 mM and 700 mM, respec- 
tively. The calculated increase in the binding energy of a 
capsid is thus 435 k B T in a surprisingly good agreement 
with experimental results. Similar results are obtained 
in model II using e.g. <7i=0.26 e/nm 2 , CT2=-0.13 e/nm 2 , 



S « 1 nm, and e p = 5 [compare Eqs. (|20"|) and ([16]) ]. Ob- 
viously, quality of the agreement crucially depends on the 
value of the surface charge density adopted, since the free 
energy scales with the square of a [see Eq.(fT5"j) and Fig. 
[3p)]. Assuming that ct=0.7 e/nm 2 , as has been estimated 
in Ref. [H (see also the data estimated in Ref. [13), the 
increase in the binding energy is about 1220 k B T, more 
than three times larger from the experimental estimate. 
Nevertheless, the order of magnitude agreement clearly 
suggest that the observed increase in the "contact en- 
ergy" between the protein subunits with the increase of 
salt concentration should be ascribed to screening of the 
repulsive electrostatic protein-protein interactions rather 
than to some salt-induced conformational change in the 
protein structure as has been assumed in Ref. 0. The 
same conclusion has been found Ref. |l5|. 



Although the asymptotic behavior of energies in the 
DH approximation [Eq. ((20] does not show the differ- 
ence in cases when the inner charge density becomes the 
outer and vice versa, numerically exact results do show a 
slight difference in the regime appropriate for viral cap- 
sids, as illustrated in Fig. [7] The free energies are smaller 
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FIG. 7: Dependence of capsid free energies when the inner 
and outer charge densities interchange. Full lines denote the 
free energies when eri = 0.5 e/nm 2 and 02 = —0.1 e/nm 2 
(<ri = —0.1 e/nm 2 and 01 = 0.5 e/nm 2 ). Dashed lines de- 
note the free energies when ai = 0.6 e/nm 2 and a 2 = —0.2 
e/nm 2 (<7i = —0.2 e/nm 2 and 02 = 0.6 e/nm 2 ). In these 
calculations, e p = 5, 8 = 1 nm, and cq = 100 mM. 



in the case when larger charge density is on the inner side 
of the capsid. Note, however, that the difference in free 
energies in the two cases is less than about 100 k B T per 
capsid, at least for the parameters chosen for the cal- 
culations displayed in Fig. \7\ This translates in about 
k B T per protein subunit, which is of the order of 10 % 
of the total energy of protein-protein contacts. Thus, 
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the inside/outside (top/bottom) asymmetry of the pro- 
tein charge may have some influence on the precision of 
the assembly so that the proper sides of proteins face 
the capsid interior. Again, preferred curvature and steric 
constraints could be more important reasons for the pre- 
cision of viral assembly. 



V. LIMITATIONS OF THE MODELS 

Our models of an empty viral capsid are approximate 
and limited in several aspects. First, we have assumed 
that the surface charges are distributed uniformly on in- 
finitely thin spheres. This is of course an oversimplifi- 
cation of the true distribution of protein charges in the 
capsid. Model II does account for the essential char- 
acteristics of the capsid charge distribution, but it still 
does not incorporate the posible polar and azimuthal in- 
homogeneities of capsid charge. In that respect, it is of 
interest to note that the R 2 behavior of the free energy 
obtained in model I for kdhR ^ 1 is obtained, at least 
on the level of the DH approximation, even when one 
accounts for the charge inhomogeneities on the surface 
of the sphere that have the icosahedral symmetry of the 
capsid 26 . Thus, it doesn't seem likely that the spatial 
distribution of protein charge could favor the preferred 
angle of protein-protein contacts. 

A second approximation in our models is the assump- 
tion of perfect sphericity of the viral capsid. The elastic- 
ity model of viral capsids by Lidmar et a&- suggests that 
larger viruses are more facetted and that their asphericity 
is larger. This is indeed a general trend that is observed 
in the experimental studied of viruses^. However, the 
strict predictions of the theory are based on its contin- 
uum limit (i.e. large T numbers) which may not hold 
for small viruses. Nevertheless, numerical experiments 
with model capsid of smaller T-number o 11 ' 38 , performed 
as detailed in Refs. [THH do show that the aspheric- 
ity in small viruses is again governed by the interplay of 
bending and stretching energies of the viral capsid and by 
the capsid radius. Elastic properties of such small cap- 
sids that are very much different from the corresponding 
properties of large viral capsids would be needed to re- 
produce a significant asphericity in small viruses, know- 
ing the asphericity of larger viruses where the continuum 
limit of the theory is expected to work fine. Thus, the as- 
sumption of small asphericity in capsids of small viruses 
is corroborated by the discrete version of the models elab- 
orated in Refs. Dilil This, however does not account 
for the corrugation of the capsid surface that is observed 
experimentally, especially in some viruses (e.g. bacte- 
riophage P4, see Ref. |22| ). This effect is related to the 
spatial distribution of the protein mass and charge, and, 
as elaborated in the previous paragraph, is neglected in 
our model. 



VI. SUMMARY 

In summary, we find that electrostatic interactions 
make important (repulsive) contribution to the energet- 
ics empty capsid binding of several hundreds of fc^T, de- 
pending strongly on the salt concentration. In highly 
salty solutions, the electrostatic interactions are effi- 
ciently screened resulting in a larger binding energy of 
the capsids, in accordance with what was experimentally 
found for Hepatitis B virus 3 -. The results of our exact nu- 
merical studies compare favorably with limiting expres- 
sions derived previously in the Debye-Hiickel approxima- 
tion and valid for sufficiently large salt concentrations^. 
The regimes in which these approximations severely fail 
have been identified. Since the electrostatic energies of 
viral capsids scale with the second power of capsid ra- 
dius, in the regime of radii typical for viruses, we find 
that simple, curvature independent expressions for the 
adhesive, attractive interaction, such as Eq. (|13p . can- 
not explain the monodispersity of self-assembled empty 
viral capsids and that certain angle-dependent interac- 
tion acting between the neighboring protein subunits is 
needed in that respect. Our study also explains that if 
the empty viral structure that is formed spontaneously 
has the same symmetry as fully infectious virus (contain- 
ing the genetic material/polyelectrolyte) then the pres- 
sure that acts on the viral capsid arises from the poly- 
electrolyte self-interaction, and the attractive interaction 
between the polyelectrolyte and the viral capsid, i.e. the 
capsid electrostatic self-repulsion is exactly counteracted 
by the curvature dependent adhesion energy. 
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APPENDIX: DETAILS OF THE SOLUTION OF 
MODEL II IN THE DEBYE-HUCKEL 
APPROXIMATION 

The potential in the interior of the capsid can in the 
DH approximation be written as 

, , . sinh(Knf/r) 

$(r)=A v - , r<R. (A.l 

r 

In the capsid material, the potential is 

$(r) = — + D, R<r<R + 6, (A.2) 
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and outside the capsid 



$(r)=B eX ^- KDHr \r>R- 



(A.3) 



Two equations for four unknown coefficients, A,B,C and 
D are obtained from the requirement of the continuity of 
potential at r = R and r = R + S. Additional two equa- 
tions are obtained by relating the discontinuity in the 
electric displacement field to the surface charge density, 



eo£ r 
dr 



<9$ 
dr 



r-R=0~ 



d<S> 

dr 



r-R=0+ 



e e r 



r -R-5 = 0- 



d<S> 
dr 



a 2 (AA) 



By solving the thus obtained equations, one obtains fairly 
complicated expressions for the unknown coefficients, 
which we simplify somewhat by introducing substitutions 

A = A/ A, B = B/A, C = C/A, D = V/A, (A.5) 

where 

A = e a e r {(e p ~ e r )K DH S 2 + e p K DH R 2 
+ [e p {l + 2n DH R)-e r {l + K DH R)]8 
+ k dh R [e r d 2 n DH + E. p R 

+ e r {l + K DH R)S}coth(K DH R)}, (A.6) 



and 



A = {e r a 1 R 2 K DH S 2 + e p [<T 1 R 2 + (j 2 {R + S) 2 ]R 
+ e r ai(l + k dh R)SR 2 } csch(K DH R), (A. 7) 



B 



{[e p (aii? 2 + a 2 (R + 6) 2 ) - e r a 2 (R + 6) 2 }S 



+ epi^R 2 + a 2 (R + S) 2 ]R 

+ e r a 2 (R + S) 2 Skdh Rcoth(n dhR)} 
x exp[K DH (R + 5)] , 



(A. 



C = e r {cn5 2 K DH R 2 + [axR? + a 2 (R + S) 2 ]R 
+ <t 1 R 2 {1 + 2k dh R)5 

- <r 2 (R + S) 2 K DH coth(n DH R)} , (A.9) 



V = e^R 2 + a 2 (R + S) 2 } 

- e r [dii? 2 + a 2 (R + 5) 2 + a 1 R 2 K DH {R + S)] 
+ e r K DH a 2 (R + 5) 2 Rcoth(K DH R). (A. 10) 

When <5 = 0, the potential reduces to that of a single 
shell of charge with the effective surface charge density of 
er = o"i +(72 [in the DH approximation, see Eqs. (| 10[) and 
([IT])]. When kdhR 3> 1, <$ <C -ft, and e r > e p , which is a 
regime of interest for viral capsids (denoted by subscript 
vir in the equations below), the following limiting forms 
for A, A, and B coefficients apply: 

lim A = e e r [nDHR 2 {^ r HDH5 + 2e p )] , (A. 11) 



lim.4 = R 3 [(cti + a 2 )e p + aie r K DH 5} csc1i(kdhR.), 

vir 

(A.12) 

limS = R 3 [(a i + a 2 )e p + (J2€rKDH$] exp(KDifR)- 

vir 

These three limiting forms are sufficient to derive Eq. 
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